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ABSTRACT 

Apollonian circle packings arise by repeatedly filling the interstices between four mutually tan- 
gent circles with further tangent circles. We observe that there exist Apollonian packings which 
have strong integrality properties, in which all circles in the packing have integer curvatures 
and rational centers such that (curvature) x (center) is an integer vector. This series of papers 
explain such properties. 

A Descartes configuration is a set of four mutually tangent circles with disjoint interiors. An 
Apollonian circle packing can be described in terms of the Descartes configuration it contains. 
We describe the space of all ordered, oriented Descartes configurations using a coordinate 
system A^p consisting of those 4x4 real matrices W with W T Q/}W = Q\y where Qd is the 
matrix of the Descartes quadratic form Qrj = x\ + x\ + x\ + x\ — \{x\ + X2 + x% + X4) 2 and 
Qiy of the quadratic form Qw = —8x1X2 + 2x| + 2x\. On the parameter space M.jd the group 
Aut(Qf)) acts on the left, and Aut(Qw) acts on the right, giving two different "geometric" 
actions. Both these groups are isomorphic to the Lorentz group 0(3,1). The right action 
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of Aut(Qw) (essentially) corresponds to Mobius transformations acting on the underlying 
Euclidean space M 2 while the left action of Aut(Qo) is defined only on the parameter space. 

We observe that the Descartes configurations in each Apollonian packing form an orbit 
of a single Descartes configuration under a certain finitely generated discrete subgroup of 
Aut(Qr>), which we call the Apollonian group. This group consists of 4 x 4 integer matrices, 
and its integrality properties lead to the integrality properties observed in some Apollonian 
circle packings. 

We introduce two more related finitely generated groups in Aut(Qo), the dual Apollo- 
nian group produced from the Apollonian group by a "duality" conjugation, and the super- 
Apollonian group which is the group generated by the Apollonian and dual Apollonian groups 
together. These groups also consist of integer 4x4 matrices. We show these groups are 
hyperbolic Coxeter groups. 

Keywords: Circle packings, Apollonian circles, Diophantine equations, Lorentz group, Coxeter 
group 
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1. Introduction 

An Apollonian circle packing is a packing of circles arising by repeatedly filling the interstices 
between four mutually tangent circles with further tangent circles. We call an initial arrange- 
ment of four mutually tangent circles with distinct tangents (necessarily 6 of them) a Descartes 
configuration. 




Figure 1: An Apollonian packing, labelled with circle curvatures 

Starting from any Descartes configuration, we can recursively construct an infinite circle 
packing of the Euclidean plane, in which new circles are added which are tangent to three of 
the circles that have already been placed and have interiors disjoint from any of them. The 
infinite packing obtained in the limit of adding all possible such circles is called an Apollonian 
packing. The new circles added at each stage can be obtained using Mobius transformations 
of Descartes configurations in the partial packing. 

An Apollonian packing is pictured in Figure Q in which each circle is labelled by its 
curvature, which is the inverse of its radius. The initial configuration consists of two circles of 
radius ^ inscribed in a circle of radius 1, the latter being assigned negative curvature —1 by 
a convention given in §3. This particular Apollonian packing has the special property that all 
circles in the packing have integer curvatures. We call a packing with this property integral. 
More remarkably, if one regards the circle centers as complex numbers, then one can place the 
initial circles so that every circle in the packing has "curvature x center" a Gaussian integer 
(element of the ring Z[i].) This occurs for example when the center of the outer circle is placed 



1 



at z = and the two circles of radius g have centers at z = — ^ and at z = 5. We call a 
packing with this extra property strongly integral. 

The object of this paper is to give a geometric explanation of the origin of these integrality 
properties involving both the curvatures and the circle centers. This is based on five facts: 

(1) An Apollonian packing can be described in terms of the Descartes configurations it 
contains. 

(2) There is a coordinate representation of the set A4b of all ordered, oriented Descartes 
configurations as a six-dimensional real-algebraic variety. It consists of the set of 4 x 4 real 
matrices W satisfying a system of quadratic equations, W^Q^W = Qw, which state that 
W conjugates the Descartes quadratic form Qd = x\ + x\ + x\ + x\ — \{x\ + £2 + X3 4- X4) 2 
to a quadratic form Qw = —8x1x2 + 2x| + 2x\ that we call the Wilker quadratic form, after 
Wilker )42| . This coordinate system we call "augmented curvature-center coordinates", as it 
encodes the curvatures and centers of the circles in the configuration. We term A4o with these 
coordinates the parameter space of ordered, oriented Descartes configurations. 

(3) The variety Mj$ is a principal homogeneous space for the Lorentz group 0(3, 1), under 
both a left and a right action of this group, realized as Aut(Qo) and Aut(Qw), respectively. 
The right action corresponds to Mobius transformations acting on the plane, while the left 
action acts only on the Descartes configuration space. 

(4) There is a discrete subgroup of the left action by Aut(Qu), the Apollonian group A, 
having the property that the (unordered, unoriented) Descartes configurations in any Apollo- 
nian packing are described by a single orbit of this group. If the Descartes configurations are 
regarded as ordered and oriented, then exactly 48 such orbits correspond to each packing, each 
one containing a copy of each unordered, unoriented Descartes configuration in the packing. 

(5) The Apollonian group A consists of integer 4x4 matrices. 

The last property (5) explains the existence of Apollonian packings having integrality prop- 
erties of both curvatures and centers. If an initial Descartes configuration has integral aug- 
mented curvature-center coordinates, then the same property persists for all Descartes config- 
urations in the packing, hence for all circles in the packing. The packing pictured in Figure ^ 
has this property. 

The observation that there are Apollonian packings having all curvatures integral is an old 
one, and was noted by F. Soddy [HZ], |3HJ - who also gave an extension to three dimensions, 
the "bowl of integers" . The existence of an integral matrix group explaining such curvatures 
was studied in Soderberg [22] in 1992. There are many other mentions in the literature. What 
we add in this direction is the observation that the circle centers can simultaneously have 
integrality properties. 

A main contribution of these papers is the viewpoint that these properties should be under- 
stood in terms of group actions on the parameter space Aio of all ordered, oriented Descartes 
configurations. Besides introducing coordinate systems for this parameter space (given first in 
[2"H] by three of the authors), we study the relevant discrete group actions in detail. In partic- 
ular we introduce a larger group of integer matrices acting on the left, the super-Apollonian 
group A s , whose generators have natural geometric interpretations in terms of their action on 
Descartes configurations. We prove that this group is a hyperbolic Coxeter group. 

The detailed contents of part I are summarized in the next section. Below we briefly 
indicate the contents of parts II and III. 

In part II we study the integrality properties of the Apollonian packing in more detail, and 
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their relation to the super-Apollonian group. We show that every integral Apollonian packing 
can be transformed to a strongly integral one by a Euclidean motion. We introduce super- 
packings as orbits of the super- Apollonian group, starting from a given Descartes configuration, 
and geometric super-packings which consist of the set of circles in all these Descartes config- 
urations. We show that there are exactly 8 different primitive 2 strongly integral geometric 
super-packings, and that each of these contains a copy of every primitive integral Apollonian 
packing. We characterize the set of all (primitive) strongly integral Descartes configurations 
as a collection of 384 orbits of the super-Apollonian group. 

In part III we consider to what extent the results proved in parts I and II extend to higher 
dimensions. In all dimensions n > 3 there are analogues of the parameter space of Descartes 
configurations, of left and right group actions by 0(n + 1,1), and of the Apollonian, dual 
Apollonian and super-Apollonian groups. However in dimensions n > 4 the Apollonian group 
action does not correspond to a sphere packing; the spheres will overlap. Nevertheless one can 
still study the orbits of these group actions on the parameter space. The Apollonian group in 
higher dimensions has rational entries, rather than integer entries. We show that configurations 
having all curvatures rational exist in all dimensions, and having curvature x centers rational 
exist in certain dimensions only. 

The general framework of these papers was developed by the second author (JCL), who 
also did much of the writing. This paper is an extensively revised version of a preprint written 
in 2000, which adds some new results in §4 and S5. 

Acknowledgments. The authors thank Ludwig Balke for suggesting that the two group 
actions in §3 should correspond to a left action and a right action. The authors are grateful for 
helpful comments from Andrew Odlyzko, Eric Rains, Jim Reeds and Neil Sloane. The authors 
thank the reviewer for incisive comments leading to reorganization of the paper. 

2. Summary of Results 

There are at least three ways to describe the Apollonian packing V containing a given Descartes 
configuration V. 

(Gl) [Geometric] An Apollonian packing Vv is a set of circles in C = M 2 U {oo}, which 
comprise four orbits under the action of a discrete group G^(T>) of Mobius transformations 
inside the conformal group M6b(2). The discrete group G_a(V) depends on V. 

(G2) [Algebraic] The collection of all (ordered, oriented) Descartes configurations in the 
packing Vt> form 48 orbits of a discrete group A, the Apollonian group, contained in the group 
Aut(Aio) = Aut(Qo) of left-automorphisms of the parameter space A4b of Descartes config- 
urations. The discrete group A is independent of the configuration V. 

(G3) [Holographic] The open disks comprising the interiors of the circles in the packing are 
the complement C\Ax> of the limit set Ax> of a certain Schottky group Sjy acting on hyperbolic 

2 An integral packing is primitive if the greatest common divisor of all curvatures of circles in the packing is 

1. 
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3-space H 3 , with C identified with its ideal boundary. The Schottky group Sx> depends on V. 

In this paper we will mainly consider viewpoints (Gl) and (G2). Viewpoint (G3) is de- 
scribed in Mumford, Series and Wright |33[ Chapter 7], and we treat it in Appendix B. It 
is termed "holographic" because it views the limit set of the packing as the boundary of a 
higher-dimensional object, which in principle gives information about it. In Appendix B we 
point out a connection of viewpoint (G3) to viewpoint (Gl). The term "algebraic" for (G2) 
refers to the group action being on a real-algebraic variety of Descartes configurations. 

The main emphasis of this series of papers is to study Apollonian packings in terms of 
the Descartes configurations they contain. That is, we study the packing as a collection of 
points inside the parameter space of all ordered, oriented Descartes configurations. There 
are two different group actions on this parameter space, which are a right action associated 
to viewpoint (Gl) and a left action corresponding to viewpoint (G2), with the groups both 
isomorphic to the Lorentz group 0(3, 1), a real Lie group, as explained in §3. In particular 
group actions are independent and mutually commute. The discrete group Gj^(T>) above is 
contained in the right action and the discrete group A, the Apollonian group, is contained in 
the left action. Thus viewpoints (Gl) and (G2) are complementary and coexist simultaneously 
on the space of all Descartes configurations. 

As stated in the introduction, the integer structures in the curvatures and centers of some 
Apollonian packings can be explained in terms of the viewpoint (G2), using the discrete group 
A, the Apollonian group. This viewpoint can be traced back to the "inversive crystal" in 
Wilker 021 §14]. 

In §3 we coordinatize the space M.o of (ordered, oriented) Descartes configurations as in 
|28j and describe two group actions on this space. In §3.1 we give two coordinate systems. 
The first labels such a Descartes configuration with a 4 x 3 matrix Mp, called curvature- 
center coordinates, and the second with a 4 x 4 matrix Wp, called augmented curvature- center 
coordinates. These coordinates are characterized by quadratic relations which generalize the 
Descartes circle theorem. Theorem 13.11 exactly characterizes these relations for curvature- 
center coordinates, strengthening Theorem 3.2 in |28| by formulating and proving a converse. 
For the augmented curvature-center coordinates these quadratic relations (proved in 28 ) take 
the form 

WlQ D Wv = Qw, 

which gives a conjugacy of the Descartes quadratic form Qjj to the Wilker quadratic form Qw 
These forms are indefinite of signature (3,1) and their (real) automorphism groups Aut(Qr>) 
and Aut(Qw) are isomorphic to the Lorentz group 0(3, 1). 

In §3.2 we describe a linear left action ("Lorentz action") by Aut(Qrj) and a right action 
("Mobius action") by Aut(Qy/) on the space A4o- Theorem 13.31 describes these actions. It 
shows that the space Alp is a principal homogeoneous space for the group 0(3, 1) under both 
the left action and the right action. The Mobius action is treated in more detail in Appendix A. 
In §3.3 we describe some integral elements of the Lorentz action Aut{Qo) which have geometric 
interpretations as simple transformations of a Descartes configuration. These elements are 
used in defining the Apollonian group, dual Apollonian group and super-Apollonian group 
given below. 

In §4 we describe Apollonian packings and the Apollonian group. Theorem 14 . 31 est ablishes 
the basic fact that the interiors of all circles in an Apollonian packing are disjoint. The much- 
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studied residual set A('P) of an Apollonian packing is the complement of the interiors of all 
circles; it is a set of measure zero. For later use, Theorem 15.21 gives several properties of the 
residual set A(V). Theorem 14.31 shows that the (ordered, oriented) Descartes configurations in 
a packing form 48 orbits of this group. 

In §5 we define the dual Apollonian group A 1 - and call its orbits dual Apollonian packings. 
Theorem 15.31 shows that the set of all tangency points of circles in a dual Apollonian packing 
has closure the limit set of another Apollonian packing, that generated by the dual Descartes 
configuration of any Descartes configuration generating the packing. 

In §6 we define the super- Apollonian group A to be the group generated by A and A 1 ' 
combined, and call its orbits Apollonian super-packings. It is a discrete group of integer 
matrices contained in Aut(Qjj)- Theorem 16.11 in §6.1 gives a complete presentation for A^, 
establishing that it is a hyperbolic Coxeter group. In §6.2 we add remarks on super-packings, 
which are studied at length in part II. 

In Appendix A we describe the Mobius group action in detail. Theorem 17.21 gives an 
isomorphism of this group to Aut(Qw)- 

In Appendix B we describe the Schottky group action in Mumford, Series and Wright |33j . 
We indicate some relations to the Mobius group action. 

3. Descartes Configurations and Group Actions 

3.1. Descartes Configurations and Curvature-Center Coordinates 

In 1643, Descartes found a relation between the radii for four mutually disjoint tangent circles 
of type (a) in Figure |2 below. Let ri,f2,r3,r4 be the radii of the tangent circles. Descartes 
showed 3 a result equivalent to 

1 1 1 1 1 / 1 1 1 1 \ 2 . . 

- + - + - + - = - - + - + - + - , (3.1) 

rf r% r§ r| 2 \ri r 2 r 3 r 4 / 

which is now called the Descartes circle theorem. This result can be extended to apply to all 
Descartes configurations, which can be of types (a)-(d) in the figure below. 

In Figure [2] (a) is a standard Descartes configuration in which the circles have disjoint 
interiors, and all curvatures are positive. We will also allow configurations like Figure |U(b) in 
which one circle encloses the other three. We also allow straight lines to be regarded as circles 
of curvature zero, leading to Descartes configurations like Figure [21(c) and (d), in which the 
"interior" of the "circle" defined by a straight line is defined to be a suitable half-plane. 

To extend the Descartes Circle Theorem to all Descartes configurations of types (a)-(d), we 
must define the curvatures to have appropriate signs, as follows. An oriented circle is a circle 
together with an assigned direction of unit normal vector, which can point inward or outward. 
If it has radius r then its oriented radius is r for an inward pointing normal and —r for an 
outward pointing normal. Its oriented curvature (or "signed curvature") is 1/r for an inward 
pointing normal and —1/r for an outward pointing normal. By convention, the interior of an 
oriented circle is its interior for an inward pointing normal and its exterior for an outward 
pointing normal. An oriented Descartes configuration is a Descartes configuration in which 

3 Descartes expressed his relation in a different form, obtained by clearing denominators in 13.11 . 
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(a) (b) 




(c) (d) 

Figure 2: Descartes configurations 

the orientations of the circles are compatible in the following sense: either (i) the interiors 
of all four oriented circles are disjoint, or (ii) the interiors are disjoint when all orientations 
are reversed. Each Descartes configuration has exactly two compatible orientations in this 
sense, one obtained from the other by reversing all orientations. The positive orientation of a 
Descartes configuration is the one in which the sum of the signed curvatures is positive, while 
the negative orientation is the one in which the sum of the curvatures is negative. (One can 
show that the sum of the signed curvatures is always nonzero.) Note that positive orientation 
corresponds to case (i), and negative orientation to case (ii). With these definitions, the 
Descartes circle theorem remains valid for all oriented Descartes configurations, using oriented 
curvatures. 

The geometry of Descartes configurations in Figure |2] is encoded in the curvature vector 
b = (b\, &2> ^3, bi) T where b\ is the oriented curvature of the ith circle. If Y2j=i bj > 0> then 
one of the following holds: (a) all of b%, 62,^3,^4 are positive; (b) three are positive and one 
is negative; (c) three are positive and one is zero; or (d) two are positive and equal and the 
other two are zero. 

The Descartes circle theorem gives a quadratic equation for the curvature vector b, which 
can be rewritten as 

b T Q D b = 0, (3.2) 
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in which 



Qd 



2 



1 


-1 


-1 


-1 


-1 


1 


-1 


-1 


-1 


-1 




-1 


-1 


-1 


-1 


1 



(3.3) 



is the Descartes quadratic form. Here 1 represents the vector (1, 1, 1, 1) T . 

In three of the authors of this paper showed there exist matrix extensions of the 
Descartes circle theorem which encode information about both the curvatures and the centers 
of the circles in the configuration, as follows. Given an oriented circle C with center (x, y) 
and oriented curvature b we define its curvature- center coordinates to be the 1x3 row vector 
m(C) := (b,bx,by). For the "degenerate case" of an oriented straight line H we define its 
curvature-center coordinates as 

m(H) := (0,h), 

where h = (hi, h 2 ) is the unit normal vector giving the orientation of the straight line. 

Theorem 3.1 (Extended Descartes Theorem) Given an ordered configuration T> of four 
oriented circles with oriented curvatures (b\, b 2 , b 3 , 64) and centers {(xi,yi) : 1 < % < 4}, let 
Mp be the 4x3 matrix 



bi bixi biyi 

b 2 b 2 x 2 b 2 y 2 

b3 b 3 x 3 b 3 y 3 

64 64X4 647/4 



(3.4) 



We include the "degenerate cases" where some circles are oriented lines. If this configuration 
is an oriented Descartes configuration, then M = Mp has a nonzero first column and satisfies 



M T Q D M 




2 
2 



(3.5) 



Conversely, any real solution M to 13.5\) with a nonzero first column is the curvature- center 
coordinate matrix M-p of a unique ordered, oriented Descartes configuration T>. 



Remark. The hypothesis of a nonzero first column is necessary in Theorem 13.11 Given the 
matrix M of a Descartes configuration, the matrix M obtained by zeroing out its first column 
will continue to satisfy ()3.5|) . and will not come from a Descartes configuration. 

Proof. An n-dimensional version of this theorem was stated as Theorem 3.2 in [2Sj, and the 
"if" direction of the result was proved there. (The converse part of the theorem was not proved 
there.) More precisely, in [!2H] it was proved that the curvature-center coordinate matrix M 
of any ordered, oriented Descartes configuration satisfies (|3,5|) . The first column of the matrix 
M of a Descartes configuration cannot be identically zero because at least two circles in any 
Descartes configuration have nonzero curvature. 
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To establish the converse part of the theorem, we need Theorem 13.21 below, which was 
independently proved in [2S]. We postpone the proof to the end of §3.1. ■ 

The curvature-center coordinates m(C) uniquely determine a circle in general position, but 
they fail to be unique for the degenerate case of a line, where the information (0, h) determines 
only a direction normal to the line and not its position. 

In [2S] three of the authors defined augmented curvature-center coordinates, (called ACC- 
coordinates for short), w(C) of an oriented circle C with signed curvature 6 and center (x,y), 
which resolve this ambiguity. These represent C by a row vector 

w(C) := (b,b,bxi,bx 2 ) (3.6) 

in which 6 is the signed curvature of the oriented circle obtained by inversion in the unit circle. 
The operation of inversion in the unit circle acts on M? by 

(x,y) ^ {x \y) = ( 2 X . 2 , 2 V . 2 ). 

x z + y z x z + y z 

It maps a circle C of finite (oriented) radius r to the circle C with center x = x/(|x| 2 — r 2 ) and 

- | x |2_ r .2 

oriented radius r = , J_ a , having oriented curvature b = — . If C is a straight line with 

|X| V ^ V 

specified normal direction, we determine C as the image of this line, with orientation coming 
from the specified normal. In the degenerate case that C is a straight line we define 6 = 0. In 
all cases, 

(bx,by) = * = * = (bx',by'), (3.7) 

so that C and C have the same curvature x center data. 

Augmented curvature-center coordinates provide a global coordinate system: no two dis- 
tinct oriented circles have the same coordinates. The only case to resolve is when C is a straight 
line, i.e., 6 = 0. The relation ()3.7|) shows that (6, 6x,6y) are the curvature-center coordinates 
of C, and if 6 ^ 0, this uniquely determines C; inversion in the unit circle then determines 
C. In fact, 6 in this case is twice the distance of C from the origin. In the remaining case, 
6 = 6 = and C = C is the unique line passing through the origin whose unit normal is given 
by the remaining coordinates. 

Given a collection [C\, C 2 , C3, C4) of four oriented circles (possibly lines) in M 2 , the aug- 
mented matrix W associated with it is the 4x4 matrix whose j-th row has entries given 
by the augmented curvature-center coordinates w(Cj) of the j-th circle. The following result 
characterizes oriented Descartes configurations. 

Theorem 3.2 (Augmented Euclidean Descartes Theorem) Given an ordered configu- 
ration T> of four oriented circles (or lines) {Cj : 1 < i < 4} with curvatures (b±, 62, 63, 64) 
and centers {(xi,yi) : 1 < i < 4}, let Wx> be the 4x4 matrix 



W v 



61 61 61x1 hyi 

6 2 6 2 6 2 x 2 6 2 i/2 

63 63 63x3 b 3 y 3 

64 64 64X4 643/4 



(3.8) 



S 



IfV is an (ordered, oriented) Descartes configuration then\Vx> satisfies 



W T Q D W 






-4 








-4 

















2 














2 



(3.9) 



Conversely, any real solution W to \3. iJ\) is the augmented matrix W© of a unique ordered, 
oriented Descartes configuration T>. 

Proof. This is proved as the two-dimensional case of Theorem 3.3 in [2H]. ■ 
We call the quadratic form Q\y defined by the matrix 



Qw = 





-4 





-4 



2 

2 



(3.10) 



the Wilker quadratic form. We name this quadratic form after J. B. Wilker [12], who in- 
troduced in spherical geometry a coordinate system analogous to augmented curvature-center 
coordinates, see [321 §2 p. 388-390 and §9]. However Wilker did not explicitly formulate any 
result exhibiting the quadratic form Qw, see [23 p. 349 ]. 

Theorem 13.21 identifies the set of all ordered, oriented Descartes configurations D with the 
set A4o of real solutions W = to the matrix equation (|3.9|) , This equation states that 
the augmented matrix coordinates of an oriented Descartes configuration give an intertwining 
map between the Descartes form and the Wilker form. The set Mo has the structure of a 
six-dimensional affine real-algebraic variety. 

It is possible to refine the parameter space to a moduli space Aio of (unordered, unoriented) 
Descartes configurations as an orbifold Aio/ ~ obtained by quotienting by a finite group of 
order 48 (acting on the left). This group is generated by the 4x4 permutation matrices 
(permuting rows) and —I, which reverses total orientation. This orbifold has singular points 
at those ordered, oriented Descartes configurations which remain invariant under a non-trivial 
permutation matrix. For our purposes it is more convenient to use the parameter space A4j$ 
which is a smooth manifold. 

Both the Descartes quadratic form and the Wilker quadratic form are equivalent over the 
real numbers to the Lorentz quadratic form Qc(x) :- 
matrix 



x o + x i + x 2 + x h w hh associated 



Qc 



-10 
10 
10 
1 



(3.11) 



For any quadratic form Q, let Aut{Q) be the group of automorphisms under the "congruence 
action" , defined by 

Aut(Q) = {U g GL(4,R) : U T QU = Q}, 
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where Q is the symmetric matrix that represents Q. The Lorentz quadratic form has a large 
group of automorphisms Aut(Qc), which is exactly the real Lorentz group 0(3, 1). 

The Descartes form and Wilker form are not only equivalent to the Lorentz form over the 
real numbers, but also over the rational numbers. For the Descartes form one has 



Qd = Jo Q-c-Jo, 



in which 



J = - 

u 2 



1 1 
1 -1 
1 -1 
1 1 



(3.12) 



(3.13) 



and Jn = Jn = Jn . It follows that 



Aut(Q D )=3o 1 0(3,l)3 . 



(3.14) 



The rational equivalence of the Wilker quadratic form Qw to the Lorentz form Qc fol- 
lows from Theorem 13.21 as soon as we exhibit a Descartes configuration T> whose augmented 
curvature-center coordinates Wx> are a rational matrix. One is given by 



W = Wd 



2 

2 

1 

1 



1 

-1 

1 
-1 



(3.15) 



It corresponds to the positively oriented Descartes configuration pictured in Figure |31 in which 
the dotted line is the x-axis and the two circles touch at the origin. 







1 



1 







Figure 3: A strongly integral Descartes configuration 



We then have 



Q w = A T Q £ A, 



(3.16) 
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with 



JqWq 



1 

-1 









1 1 
-1 1 



(3.17) 



It follows that 

Aut{Q w ) = A _1 0(3,1)A. 
Finally we note that the determinants of these quadratic forms are given by 
det(Q D ) = det(Q £ ) = -1 and det(Q W /) = -64. 



(3.18) 



(3.19) 



Taking determinants in ()3.9|) yields that augmented curvature-center coordinate matrices W-p 
have 



det(Wx») = ±8. 
Examples show that both values ±8 occur. 



(3.20) 



Proof of the "only if" part of Theorem 13.11 We note that if M is the curvature-center 
coordinate matrix of a Descartes configuration, this configuration is unique. Indeed at least 
two circles in the configuration have nonzero curvature, so are determined by their curvature- 
center coordinates. This gives enough information to find the location of any lines in the 
configuration, and their signed normal vectors are determined uniquely by their curvature- 
center coordinates. Thus it suffices to show that there exists a 4 x 4 matrix W whose last 
three columns agree with M, such that 



W T Q D W = Q w 






-4 








-4 

















2 














2 



(3.21) 



If so, then Theorem 13 . 21 implies that W is the augmented curvature-center coordinate matrix of 
a unique ordered, oriented Descartes configuration T>. It follows that M is the curvature-center 
coordinate matrix of T>. 

This discussion implies that the extension W is unique, and we proceed to find it by 
transforming the problem from the Descartes form Qd to the Lorentz form Qc = JqQdJq. 
This proof seems "pulled out of thin air" in that it guesses to do the "translation" by the matrix 
T below, in terms of which the extended matrix W has a simple form. We set M = JoM with 
Jo given by (|3.13|) . so that Jq = Jq = Jq 1 - Then 



M T Q £ M 



M T Q D M 



Qo = 




2 
2 



(3.22) 



We now note that 



Mn = ^( m n + m i2 + mi3 + mu) / 0. 
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To verify this, note that the top left entry of (|3.5|) gives 

Qd (ran, mi2, mi3, mu) = (mil + m? 2 + m? 3 + mf 4 ) - - (mu + m i2 + m i3 + m u ) 2 = 0, 

and if mu + m\2 + 77113 + mi4 = then we would obtain mfj + mf 2 + + m^ 4 = 0, which 
forces all m\j = 0, contradicting the hypothesis. Now there exists a unique matrix 

a b c 
1 
1 

such that the first row of MT is [1,0,0], which takes a = (Mu)" 1 , so that T is invertible. We 
write 

1 

V 2 1 V 2 2\f2 W 2 3\/2 
«31 V 32 V2 V 33 \[2 
«41 V42V2 V4 3 V2 



MT 



and then have 



T T M T Q £ MT = T T Q T = Q 




2 
2 



(3.23) 



This matrix equation is equivalent to the assertion that 



V 



V21 V22 V23 
V31 v 32 v 33 

7J41 V42 ^43 



is an orthogonal matrix, i.e. V T V = VV T = I. We now define 



W 



2 10 

-2v 2 i V21 V22V2 V23V2 

-2v 31 u 3 i v 32 V2 v 33 \/2 

-2V4i U41 V42V2 



whose last three columns match MT, and verify by direct calculation using (|3.23j) that 

W T Q £ W = Q w . 
We define W := J^W = J W, and find that 

W T Q D W = W T J Q£J W = Q w . 
This is the desired lift, since the last three columns of W are exactly M. ■ 
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3.2. Mobius and Lorentz group actions 

The augmented Euclidean Descartes theorem immediately yields two group actions on the 
space Aio of ordered, oriented Descartes configurations. The group Aut{Qo) acts on the left 
and the group Aut{Qw) acts on the right, as 

Wp^UW^V" 1 , with U G Aut(Qo), V G Aut(Qw)- 

The two group actions clearly commute with each other. Both groups are conjugate to the 
real Lorentz group 0(3, 1), and therefore they each have four connected components. These 
components are specified for any Y G 0(3, 1) by the sign of det(Y), and by the sign of "total 
orientation" , which is the sign of Yn = efYei, in which ef = (1,0,0,0). For U G Aut(Q D ), 
it is the sign of 1 T U1, in which 1 T = (1,1,1,1); and for V G Aut(Qw), it is the sign of 
ei in which A is any matrix satisfying 0(3, 1) = AAut(Qw)A. , such as (|3.17j) . 
The parameter space A4o also has four connected components, specified by similar invari- 
ants, which are the sign of det(W) and the total orientation, which is the sign of the sum 4 of 
the (signed) curvatures of the four circles in the Descartes configuration. We let 

Mo = M\ U ML U m\_ U M^l, 

in which the subscript describes the sign of the determinant and the superscript the orientation, 
with | being positive orientation. We also let 

denote the set of positively oriented Descartes configurations. The relevance of this decompo- 
sition is that the Apollonian group defined in §4 leaves the set of positively oriented Descartes 
configurations ML invariant. 

The action on the right by the elements of Aut{Qw) maps circles to circles, since the 
circles in a Descartes configuration correspond to the rows in the matrix Wx> of an ordered, 
oriented Descartes configuration. This right action can essentially be identified with the gen- 
eral Mobius group M6b(2) generated by the linear fractional transformations acting on the 
one-point compactification R 2 of R 2 , which is PSL(2,C) = SL(2, C)/{±I}, together with 
complex conjugation z z. and has two connected components. More precisely, Aut(Q\y) is 
isomorphic to the direct sum of the Mobius group with a group of order two, which has four 
connected components, as explained in Appendix A. 

The action on the left, by Aut(QD), produces a new oriented Descartes configuration 
whose coordinates mix together the coordinates of the different circles in the original Descartes 
configuration. This group action does not make sense as an action on individual circles in 
the configuration. This group action is intrinsically associated to the 6-dimensional space of 
oriented Descartes configurations. 

Theorem 3.3. (1) The groups Aut{Qo) and Aut{Qw) «re conjugate to Aut{Qc) = 0(3, 1). 

(2) The group Aut^Qu) acts transitively on the left on the space A4b of all ordered, oriented 
Descartes configurations. Given two such Descartes configurations T> and T> there exists a 
unique U G Aut{Qo) such that UWx> = W-p'. 

4 The sum of the signed curvatures of the circles of an oriented Descartes configuration cannot be zero. 
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(3) The group Aut{Qw) acts transitively on the right on the space of all ordered, oriented 
Descartes configurations Mo- Given two such Descartes configurations T> and T> there exists 
a unique V G Aut{Q w ) such that WpV" 1 = W v >. 

(4) The action of Aut(Qo) on the space Mo commutes with the action of Aut(Qw)- 

Remark. The left action by Aut(Qo) and the right action by Aut{Qw) on Mo given by The- 
orem I3~3"l can be identified with a left and right action of the Lorentz group 0(3, 1) = Aut(Qc) 
on Aio, using 1)3. 14j) and (|3.18j) . respectively. Theorem 13.31 shows that both these actions are 
transitive, and that the stabilizer of a point is the identity element. This is equivalent to saying 
that the space Mb is a principal homogeneous space (or torsor) for 0(3, 1) for either action. 

Proof. Part (1) follows from the conjugacy between the Descartes and Wilker forms and the 
Lorentz form Qn given in (|3.12|) and (|3.16j) . respectively. These give 

Aut(Q D ) = J l Aut(Q c )J = 3^0(3, 1)J , 

and 

Aut(Q w ) = A^AutiQ^A = A _1 0(3, 1)A 

Parts (2) and (3) follow immediately from (1). Given a fixed W 
that 

Mo = Aut(Q D )W. 

Taking Wo in (|3.15|) we have 

(WW fQ fl (WW x ) = (WA- 1 J ) t Qd(WA- 1 J ) = Qd, 

so that WWq 1 G Aut(Q£>) and Aut(Q£>)WWQ 1 forms a single orbit of Aut(Q£>), The map 
Wp i — ► WjiWq 1 from A^b into Aut(Qo) is one-to-one since Wo is invertible, and it is onto 
since the domain includes Aut{Q£>yW. Thus (|3,24|) follows, and this gives (2). We similarly 
obtain 

Mo = WAutiQw)- 1 = WAut(Q w ), 

which gives (3). 

Finally, part (4) follows directly from parts (2) and (3), since in general a left and a right 
action of two matrix groups on a space of matrices commute. ■ 

In this paper we mainly study structures associated to the left action of Aut{Qo) on Mo- 
We term the left action the Lorentz action, although this is a misnomer, since the right action 
can also be identified with a Lorentz group action. However we wish to assign different names 
for the two actions, and the right action is conveniently called the Mobius action, because of 
its relation to a Mobius group action on circles on the Riemann sphere, detailed in Appendix 
A. 



Wp G Mo, we assert 
(3.24) 
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3.3. Distinguished Elements of Aut^Qo). 

We now describe some specific elements of Aut{Qo) that have a nice geometrically interpretable 
action on every Descartes configuration, visualizable in terms of inversions, and whose associ- 
ated matrices have integer (or half-integer) entries. 

The first set of four operations correspond to inversion in the circle determined by the 
three intersection points of circles in a Descartes configuration T> that avoid one particular 
circle. There are four possibilities for these. This inversion fixes the three circles involved 
in the intersections and moves the fourth circle, to the unique other circle that is tangent 
to the first three circles. For this reason we call it the reflection operator. Let S\ = S\\D\ 
denote the Mobius transformation of this kind that moves the circle C\. In particular, Si 
maps the Descartes configuration V = (Ci, C2, C3, C4) to s\ (T>) = (C[, C2, C3, C4). This 
reflection operator depends on the specific Descartes configuration. However for all Descartes 
configurations T> there holds 

W Sl(D) = SiWp, (3.25) 

where Si £ Aut(Qjj) is given by 

" -1 2 2 2 " 
10 
1 ~ 10 
1 _ 

belongs to Aut(Qo) and is independent of the location of T>. The geometric action of S\ is 
pictured in Figure 0J The other three operations 32,33,54 give similar matrices 82,83,84, 
obtained by permuting the first and j-th rows and columns of Si for 2 < j < 4. 




Figure 4: The reflection operator S\ 



To prove (|3.25fl we check that Si £ Aut(Qr>). Then for any Descartes configuration T> 
Theorem 13 . 21 gives that SiWx> = Wx>/ for some Descartes configuration T>' . The configuration 
T>' necessarily has three oriented circles fixed, with the fourth moved, so the fourth one must be 
the unique other circle tangent to the given three, and its orientation is uniquely determined 
by the other three orientations. Now W Sl (x>) corresponds to that Descartes configuration 
consisting of the same four ordered circles, which has the same (positive or negative) orientation 
as V. 

A second set of four operations corresponds to inversion in one of the four circles in a 
Descartes configuration. Let that circle be C\. Now C\ remains fixed, while the other three 
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circles change. Denoting this inversion by sj", it is easy to show that 

W sHv) = SiW v , (3.26) 

where Sf £ Aut{Qo) is given by 

-10 0" 
2 10 
2 10' 
2 1 _ 

The corresponding matrices S^", S3", S4 are obtained by permuting the first and j-th rows and 
columns of S^ for 2 < j < 4. The operation sj~ is pictured in Figure 



C 2 I C 3 J 

Figure 5: The inversion operation 

Finally we describe an operator, which we call the duality operator D, and which is an 
involution, as follows. Given a Descartes configuration T>, there exists a Descartes configuration 
T>' obtained from it which consists of the four circles each of which passes through the three 
intersection points avoiding one circle. The resulting configuration has the same six points of 
tangency as the original configuration, and the circles in the configuration T>' are perpendicular 
to the circles of T> at these tangency points. This is pictured in Figure H3 

Straightforward computation shows that 

Wzv = DW P , 

with 

-1 1 1 1 " 
1-111 
1 1-11' 
1 1 1 -1 _ 

Note that D = — e Aut(Qa)- There is a Mdbius transformation = which sends 
V i — > d(V) = V , and it depends on V. 
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Figure 6: The dual operation 



4. Apollonian Packings and the Apollonian Group 

Apollonian circle packings are infinite packings of circles recursively constructed from a given 
positively oriented Descartes configuration T>. For simplicity consider a positively oriented 
Descartes configuration T> = {C\, C%, C3, C4} in which circle C4 encloses the other three, so 
that the interior of C4 includes the point at infinity. The uncovered area consists of four 
lunes, indicated by the shading in Figure [7| 




Figure 7: A zero-th stage packing V-jy and uncovered regions (lunes) 

The zeroth stage packing vffi = T>. At the first stage we inscribe a circle in each lune, 
to obtain a circle packing containing 8 circles. Each of these circles lies in a unique 
(unordered) Descartes configuration in the first stage packing, consisting of it and the three 
circles it touches. The uncovered area that remains consists of 12 lunes. See Figure EJ 
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At the second stage we inscribe a circle in each of these lunes, which produces the second 
stage packing vffi ■ Continuing in this way, at the n-th stage we add 4 • 3 n_1 circles, and 
each of these lies in a unique Descartes configuration in the n-th stage partial packing. The 
Apollonian packing Vv associated to V is the limit packing 

oo 

Vv := (J • (4- 1 ) 

n=l 

One can regard an Apollonian packing V as a geometric object consisting of an infinite 
collection of circles. These circles are described as the four orbits of a group of Mobius trans- 
formations Ga(V) acting on the circles (C\, C2, C3, C4) in original Descartes configuration V. 
The group 

G A {V) := (8^82,83,84), (4.2) 

in which Sj is the inversion with respect to the circle that passes through the three intersection 
points of T> that do not include the circle Cj, as indicated in §3.3. One can check that the 
circles added at the n-th stage of the construction above correspond to words s^s^ • • • Sj„ of 
length n in the generators of this group, in which any two adjacent generators are distinct, i.e. 

The key geometric property of Apollonian circle packings is that they can be viewed as 
packings of disks having the circles as boundary. To each circle in the Riemann sphere cor- 
respond two disks having the circle as boundary. Recall that the "interior" of an oriented 
circle corresponds to making a choice of one of these two disks, indicated by a choice of sign of 
the curvature, when it is nonzero (and by a choice of normal vector if the curvature is zero). 
The disk packing property of an Apollonian packing corresponds to treating its Descartes 
configurations as positively oriented. 

Theorem 4.1. For any Apollonian packing generated by a positively oriented Descartes con- 
figuration, the interiors of all circles in the packing are disjoint. 
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Proof. We generate the packing from a single Descartes configuration T> chosen to have 
positive total orientation, so that the interiors of its four circles defined by the total orientation 
are disjoint. Positive total orientation is preserved by the action of Mobius transformations; 
see Theorem 17. II in Appendix A. Let 

m 
n=l 

where T> is a positively oriented Descartes configuration. We prove by induction on m the 
following statement: Vx>[m] consists of 2(3 m + 1) circles which form the boundary of 2(3 m + 1) 
disks and 4-3 m lune areas. The interiors of the disks, as well as the lunes, are mutually disjoint. 

The base case m = is easy. Since any two positively oriented Descartes configurations 
are equivalent under the Mobius transformations, which map circles to circles and preserve 
the (total) orientation, (c.f. Theorem 7.2, Appendix A), we can simply check the Descartes 
configuration given in Figure 7, which has 4 circles that form the boundary of 4 disks and 4 
lunes with disjoint interiors. 

Assume the statement holds for m > 0. We carry out the inductive step for m + 1. By the 
inductive construction above, 

Vv[m + l]=Vv[m}UVg* l) . 

The circles of Vj)[m] remain in T , T>[m + 1] with their interiors untouched. In each lune of Vj)[m] 
a new circle is inscribed, breaking the original lune into 4 parts-one disk and three smaller 
lunes, with disjoint interiors. Combining the inductive hypothesis, we have that Vx>[Tn + 1] 
consists of 2(3 m + 1) +4- 3 m = 2(3 m+1 + 1) many circles with empty interiors, and the uncovered 
area are 3 • 4 • 3 m = 4 • 3 m+1 lunes with disjoint interiors. Each new Descartes configuration 
shares the interior of three circles with the Descartes configuration generating it; this implies 
it has the same orientation as the preceding one; hence it has positive total orientation. This 
proves the statement for m + 1. ■ 



We first consider an Apollonian packing as a geometric object. A geometric Apollonian 
packing V is the point set consisting of a countable collection of circles on the Riemann sphere 
C = C U {oo} obtained by the construction above. It has Hausdorff dimension 1, since it is 
a countable union of sets of Hausdorff dimension 1 (circles). Let G(V) denote the group of 
Mobius transformations that leave V invariant. This group can be proved to be a discrete 
subgroup of the group of all Mobius transformations, which acts transitively on the circles in 
the packing. The group G(V) contains Gj^(T>) above as a subgroup of index 24, with cosets 
given by 24 Mobius transformations whose effect is to fix a generating Descartes configuration 
T> and permute the four circles in it. (These are nontrivial facts, and we do not prove them 
here.) 

The geometric Apollonian packing is not a closed set on the Riemann sphere. We define 
A(V) to be its closure on the Riemann sphere, and call it the residual set of the Apollonian 
packing. These sets are prototypical examples of fractal sets, and have been much studied; in 
1967 Hirst [2£] showed these sets have Hausdorff dimension strictly between 1 and 2. For later 
use, we summarize properties of A(V) in the following result. 
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Theorem 4.2. The residual set A(V) of a geometric Apollonian packing V has the following 
properties. 

(1) A(V) is the complement in the Riemann sphere of the interiors of all circles in the 
packing. Here "interior" is defined by a positive orientation of a Descartes configuration in 
the packing. 

(2) A(V) has a Hausdorff dimension which is independent of the packing V , that satisfies 
the bounds 

1.300197 < dim H (A(T)) < 1.314534. 

In particular, A(V) has Lebesgue measure zero. 

(3) A(V) is the closure of the countable set of all tangency points of circles in the packing. 

(4) A(V) is invariant under the action of the discrete group G(V) of Mobius transforma- 
tions. 

Remark. Property (1) is the more usual definition of residual set of an Apollonian packing. 

Proof. Let I(V) be the open set which is the union of the interiors of all circles in the 
packing. It is clear that 

A(V) C A*(V) := C\I(V). 

The set A* (V) is the usual definition of the residual set of an Apollonian packing; we will show 
A(V) =A*(V). 

The result of Hurst |25j that the Hausdorff dimension of A* (V) is strictly less than 2 implies 
that A*(7 3 ) has Lebesgue measure zero. The sharper bounds on the Hausdorff dimension 
stated here were obtained in 1973 by Boyd [Tj- A later result of Boyd [5] suggests that 
dim//(A*("P)) = 1.3056 db .0001. Thomas and DharjlHl give a non-rigorous approximation 
scheme suggesting that this Hausdorff dimension is 1.30568673 with an error of 1 in the last 
digit. See Falconer [TBI pp. 125-131] for proofs of 1< dim#(A*(P)) < 1.432. 

Now let z G A*(V). The interiors of the circles in the Apollonian packing cover all but a 
measure zero area in a neighborhood of z. This requires infinitely many circles, so their radii 
must go to zero, since there are only finitely many circles of radius larger than any positive 
bound. Since each of these circles has a tangency point on it, the point z is a limit point 
of such tangency points. Since the tangency points are contained in V, so we conclude that 
A*(V) C A{V). This gives A(V) = A*(V), at which point properties (1), (2) and (3) follow. 

Property (4) follows by observing that the group invariance of V carries over to its closure, 
by applying it to any Cauchy sequence of points in V . ■. 

The main viewpoint of this paper is to treat an Apollonian circle packing V as described 
"algebraically" by the set of all Descartes configurations it contains, which we denote In 
the construction process above, except for the four circles in T>, each other circle C in the pack- 
ing corresponds to the unique (unordered) Descartes configuration containing it which occurs 
at the stage of the construction where C first appears. Each (unordered, unoriented) Descartes 
configuration appears exactly once in this correspondence, except for the base configuration T> 
which corresponds to the initial four circles. These Descartes configurations can be described 
in terms of orbits of a discrete subgroup A of Aut(Qx>) called here the Apollonian group. 
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The Apollonian group A is a subgroup of Aut(Qo) denned by the action of the reflection 
operations on Descartes configurations described in §3.3. These correspond to inversions in 
the circles passing through the three tangency points in a Descartes configuration that do not 
include one fixed circle in the configuration. 



Definition 4.1. The Apollonian group A is the subgroup of Aut(Qo) defined by 

A := (Si, S2, S3, S4) , 



(4.3) 



where 
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(4.4) 



We now characterize the set of ordered, oriented Descartes configurations fS>(V) in an Apol- 
lonian circle packing V in terms of this group. We note that a single unordered, unoriented 
Descartes configuration, corresponds to exactly 48 ordered, oriented Descartes configuration, 
since there are 24 choices of ordering of the four circles, and two choices of orientation. 

Theorem 4.3. The ordered, oriented Descartes configurations in an Apollonian circle packing 
comprise a union of 48 orbits of the Apollonian group. Each of these orbits contains exactly 
one ordered, oriented representative of each (unordered, unoriented) Descartes configuration 
in the packing. 



Proof. Let V = (C±, C2, C3, C4) be the zeroth stage packing of V, which is ordered and 
positively oriented. It is enough to show that A(T>) contains exactly one ordered, positively 
oriented representative of each unordered, unoriented Descartes configuration in the packing. 
Note that Sf = S| = S3 = S4 = I, and for any word S^S^ . . . Sj„, (ij 7^ ij+i), the Descartes 
configuration V = (C[, C' 2 , C 3 , C' 4 ) with ACC coordinates Wp< = Sj 1 Sj 2 . . . Sj n Wx> is a posi- 
tively oriented Descartes configuration containing a unique circle C[ in the n-th stage partial 
packing. Since no two circles at level n touch for n > 1, by induction it is easy to see that 

A = {S h S i2 . . . S in I n > Q,ij / ij+i}, 

and 

A{V) = {V I W v , = S h S i2 . . . S in W v , n > 0,ij + i j+1 }. 

Furthermore, for any ordered, positively oriented configuration V" in V, there is a unique 
permutation a such that o~(T>") = D' £ A(T>), where there is a unique sequence i±i2 . . . ,i n , 
1 < ij < 4 and ij 7^ ij+i, such that Wp/ = S^S^ . . . Sj n Wx>. This proves that all the ordered, 
positively oriented Descartes configurations in an Apollonian circle packing comprise a union 
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of 24 orbits of the Apollonian group. The theorem follows by counting both the positively and 
negatively oriented Descartes configurations. ■ 

Theorem 14 . 31 shows that the orbit under the Apollonian group of a single (ordered, oriented) 
Descartes configuration completely describes an Apollonian packing. 



5. Dual Apollonian Group 



As explained in §3, the operation of inversion in the each of the individual circles of a Descartes 
configuration is described by an integral matrix in Aut(Qjj). This leads us to the following 
definition. 

Definition 5.1. The dual Apollonian group A 1 " is the subgroup of Aut(Qr>) generated by the 
matrices 
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(5.5) 



We note that = Sj , the transpose o/Sj. 

By analogy with Theorem 14.31 we might think of an orbit of the dual Apollonian group 
acting on a single oriented Descartes configuration D, as a "dual Apollonian circle packing". It 
can be viewed algebraically as a set of Descartes configurations or geometrically as a collection 
of circles in the plane. In the algebraic viewpoint the orbit is a discrete set of points in the 
parameter space Mj$ so is a discrete object. Prom the geometrical viewpoint, as a collection 
of circles, which we denote V^, it has the following weak "packing" property. 

Theorem 5.1. No two circles in a dual Apollonian circle packing cross each other. That is, 
two circles in distinct Descartes configurations in a dual Apollonian circle packing either 
coincide, or are tangent to each other, or are disjoint. 

We do not give a proof of this theorem here, as it follows from a similar result proved for 
the super-Apollonian packing in §3 of part II. It also is a special case of an n-dimensional 
generalization proved in §4.2 of part III. 

Figure ED pictures the circles in a dual packing V 1 " . These are circles in Descartes con- 
figurations dual to those in the Apollonian packing in Figure ^ The fractal-like part of this 
figure is the closure of the set of points at which circles in the dual packing touch. It can be 
shown that infinitely many circles touch at each such tangency point. Below we shall show 
that this fractal set coincides with the limit set of the Apollonian packing generated by the 
dual Descartes configuration T)(Vq). 

The dual Apollonian group has a simple relation to the Apollonian group. 
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Figure 9: A dual Apollonian packing 
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Theorem 5.2. The dual Apollonian group A is conjugate to the Apollonian group A using 
the duality operator D 6 Aut{Qo). This holds at the level of generators, with 

D T SiD = for 1 < i < 4. 

where D = D 1 = D T . 

Proof. This is a straightforward computation. ■ 

Theorem 5.3. Given a dual Apollonian packing generated by Descartes configuration 
T>q, let £ be the set of all intersection points of the circles in A ± (T>q). Then the closure A := £ 
is equal to the residual set A(7 : 'd(Do)) °f ^ e Apollonian packing generated by the dual 

Descartes configuration D(2?o)- 

Proof. Theorem I4.2f 3) states that the residual set A('P) of a geometric Apollonian packing 
V is the closure of the set of tangency points £' in V ■ Thus it will suffice to show that 
the set of tangency points £ of the circles in the dual packing V^ Q coincides with the set of 
tangency points £' of circles in the Apollonian packing Pd(z>o)- Each intersection point of two 
circles in the dual Apollonian packing V^ Q is an intersection point of circles in some Descartes 
configuration T> of A ± (T>q). (This follows from the recursive construction of the dual packing.) 
Each intersection point or circles in T> is an intersection point of circles in the dual Descartes 
configuration D(2?), which is the configuration having augmented curvature-center coordinates 
DWx>. Now D(D) belongs to the Apollonian packing A(D(Dq)) because 

DW D = D(S,^ • • • Si m W Vo ) = S n • • • S im (DW Vo ), 

using Theorem 15.21 Thus the set £ of intersection points of circles in the dual Apollonian 
packing is contained in the set £' of intersection points of circles in the Apollonian packing 
Vt)(p \. The argument reverses to show the converse is also true: each intersection point of 
circles in the Apollonian packing "Pd(x> ) ^ s an intersection point of circles in the dual Apollonian 
packing V^ Q . Thus S = ■ 

We have the following dichotomy between the geometric and algebraic views of Apollonian 
packings and dual Apollonian packings. Viewed geometrically as collections of circles, a dual 
Apollonian packing V 1 - appears quite different from an Apollonian packing, as the circles are 
nested to an infinite depth. However when viewed algebraically as a collection of ordered, 
oriented Descartes configurations B(7 3_L ), it is a discrete set given by the orbit of a group A 
conjugate to the Apollonian group S, and so is a similar object. 

6. Super- Apollonian group 

The Apollonian group A and dual Apollonian group A^ both consist of integer automorphs 
in Aut(Qr))- We obtain a larger integral group by combining the two groups, as follows. 

Definition 6.1. The super- Apollonian group A s is the subgroup of Aut(Qrj) generated by 
the Apollonian group A and A 1 - together. We have 

A s = (Si, S2, S3, S4, Si, S2 , S^, Sf ). 
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The super- Apollonian group A s is a discrete subgroup of Aut(Qr>) as a consequence of the 
fact that all its members are integral matrices. 

6.1. Presentation of super- Apollonian group 

We now determine a presentation of the super- Apollonian group, showing that it is a hyperbolic 
Coxeter group. This result implies that both the Apollonian group and dual Apollonian group 
are finitely presented, and are hyperbolic Coxeter groups. 

Theorem 6.1. The super- Apollonian group 

A s = (Si, S2, S3, S4, Sf, S2 , S3, Sf ) 

is a hyperbolic Coxeter group whose complete set of Coxeter relations are 

S? = (S^) 2 =1 for 1 < 2 < 4. (6.1) 

(S^) 2 = (S+Si) 2 = I if i^j. (6.2) 

Proof. Recall that Sf- = Sf. The group A s satisfies the relations S? = (Sf ) 2 = I for 
1 < i < 4 and it satisfies 

S«Sj = SjS< for i^j , (6.3) 

which is equivalent to ()6.2|) . 

Write words U in A s as U = U n U n -i • • • Ui, in which each Uj is a generator Sj or Sj, 
with the empty word (n = 0) being the identity element. We do not need inverses since each 
generator is its own inverse. A word U is in normal form if it has the two properties: 

(i) U fc ^ U fe _i for 2 < k < n. 

(ii) If Ufc = Sj for some j, then XJk-i 7^ Sj for all i 7^ j. 

We can reduce any word to a word in normal form, using the relations to move any symbols 
Sj as far to the right in the word as possible. We cancel any adjacent identical symbols to 
make (i) hold. Then move the rightmost Sj as far to the right as possible, as allowed by (|6.«-{j) . 
Repeat the same with the next rightmost symbol Sj. If in the process any adjacent symbols 
(Sj) 2 occur, cancel them. This process must terminate in a normal form word or the empty 
word. 

The theorem is equivalent to showing that no normal form word with n > 1 is the identity 
element in A s . This holds because the reduction to normal form used only the Coxeter 
relations, so if A s satisfies an additional nontrivial relation, there would exist some nontrivial 
normal form word that is the identity. 

We proceed by induction on the length n of a word in normal form U = U n U n _i • • • Ui. 
Note that any suffix Uj_iUj • • • Ui of U is also in normal form. For each n > 1, let X n denote 
the set of all normal form words of length n. (By a simple enumeration we can show that the 
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number of normal form words of length n is 9 • 5™ — 1. ) We measure the size of a word U 
(viewed as a 4 x 4 matrix) as 

4 4 

/(U):=l T Ul = £^l% . (6.4) 
i=l j=l 

Thus /(I) = 4. For n = 1 we have 

/(U) = /(Ui) = 8 (6.5) 
in all cases. For n > 2 if U G X n then U' = U n _iU n _2 ■ • • Ui G X n -li an d we will prove 

/(U) > /(U') . (6.6) 

If so, then 

/(U) > /(Ui) > 8, 

hence U/I, which will complete the proof. 
We let 

r(U) :=(r 1 (U),r 2 (U),r 3 (U),r 4 (U)) T = Ul, (6.7) 

be the vector of row sums of U. Now 

/(SxU') = l T S 1 r(U / ) = (-l, 3, 3, 3)r(U') = -4r 1 (U') + 3/(U / ), 
/(SfU') = l T Sf r(U') = (5, 1, 1, l)r(U)=4r 1 (U') + /(U / ), 

with similar formulas in the other cases. To prove ()6.6|) it is therefore sufficient to prove the 
following two assertions, for 1 < h < 4. 

(1) If U n = S h , so that U n _i / S h , then 2r h (U') < f(U>). 

(2) If U n = S^, so that Un-x + Sj t and J7„_i + Si for i / h, then r fc (U') > 0. 

Instead of proving (1), (2), we will prove by induction on n the following three assertions 
for all U G X n . Here (h, i, j, k) always denotes some permutation of (1, 2, 3, 4) in what follows. 

(i) For 1 < i,j < 4, with i ^ j, 

nW+TjiU) > o . 

(ii) If U n = S h , then 



r h (U) 


> 


0, 




r»(U) 


< 


r h (U) + 


rj(U)+r fc (U) 


r,(U) 


< 


r h (V) + 


r,(U)+r fc (U). 


r fe (U) 


< 


r h (V) + 


ri(U) + rj(U) 
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(iii) If U n = S£ then 

r h (U) < 0, 

n(U) < r h (U)+ rj (U)+r k (U), 

rj(U) < r h (U)+ ri (XJ)+r k (XJ), 

r fc (U) < rfcCUJ+r^UJ + r^U) . 

Note that (i) implies that at most one of the row sums of U can be negative. If proved, (i)-(iii) 
imply (1), (2), which themselves imply completing the proof of the Theorem. 

The induction hypothesis (i)-(iii) holds for n = 1, since r(Sh) is a permutation of (5, 1, 1, 1) T 
and r(S^) is a permutation of (—1,3,3, 3) T . Suppose it is true for n. For the induction step, 
write U G X n +i as U = U n+ iU' with U' G an d abbreviate r\ = rj(U'). 



Case 1. U n +\ = S/j /or some h. 

We must verify (i) and (ii) for U. We have 





= r h (S h U') 




r»(U) 


= n(S ft U') 


= r' 


r,(U) 


= ^(SftU') 




r*(U) 


= r fc (S ft U') 


= r' k 



r> h + 2r' l + 2r>+2r> k , 



To verify (i), all cases not involving r^(U) follow from the induction hypothesis. To show 
r h (U) + r;(U) > note that 

r h (U) + r;(U) = -r£ + 3r< + 2r) + 2r^ . (6.8) 

Now U n / Sh since U G X n +i, and for all cases except U n = the induction hypotheses (ii) 
and (iii) give r' h < r ■ + r'- + r' k , whence 

r ft (U) + T-i(U) > 2r^ + + = (rj + rj) + (rj + r£) > , 

using induction hypothesis (i). If U n = Sj, then by hypothesis (iii) r' h < 0, hence i J i ,r'pr' k > 
by hypothesis (i), so all terms on the right side of ()6.8|) are positive, so r^(U) + fj(U) > 0. 
The cases r^(U) + r, (U) > and r^(U) + r&(U) > follow similarly. To verify (ii) for U it 
suffices to prove 

< -r' h + 2r- + 2r'j + 2r' k , 
r'i < -r^ + 2^ + 3r^ + 3rJ., 
r j < -^A + 3r£ + 2rJ- + 3rfc, 
< -r'h + M + 3^- + 2rJ,, 



which is equivalent to 



< 2^ + 2^ + 2^, 

r h < r- + 3r^ + 37^ 

?i < 3r- + r'j + Sr' k , (6.£ 

Tfc < 3r^ + 3rj + r^. . 
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There are three cases, according as U n = Sj, U n = S/ with i ^ h and U n = Sj r In the first 
two of these, the induction hypotheses give 

r h < r i + r 'j + r L 

which with induction hypothesis (i) immediately yields (|6.9j) . If U n = then by (iii), r' h < 0, 
whence ^,r'-,rj. > by (i) so (|6.9|) is immediate. This finishes Case 1. 



Case 2. U n+ i = /or some h. 
We must verify (i) and (iii). We have 

r fc (U) = u(S£U') = -r'„ 
r,(U) = r 4 (S£U') = + 
^■(U) = r^U') = 2r' h + r' p 
r fc (U) = r fc (S£U') = 2r^ + r;, 

To prove (i), note first that 

r h (U)+r i (U)=r' h + r' i >0, 

using (i) for U', and similar inequalities hold for other cases involving r/j(U). The remaining 
cases are all of the form r,(U) + rj(U) > with i, j ^ h, and the proofs for each are similar. 
We have 

r f (U) + rj(U) = 4r' h + r • + rj- = 2r' h + (r' h + r-) + (r^ + r'j) > 2r' h , 

using induction hypothesis (i). We now show r' h > holds in all cases. Since U G X n+ \, we 
have U n 7^ and U n 7^ Sj with i ^ h. If U n = Sj for some j ^ h then hypothesis (iii) for 
U' says that r'j < 0, whence we must have r' h > 0. If U n = Sh then hypothesis (ii) for U' says 
r' h > 0. Thus r' h > in all cases, and (i) holds for U. To prove (ii) we must prove 

r' h <0, 

2r' h + r' i <3r' h + r' j + r' k , 
2r' h + r' j <3r' h + r'+r> k , 

K + r k<K + T i + r j> 

which is equivalent to 

r' h >0, 

<<r' h + r'+r' k , 

r 'j < r 'h + < + r 'k> 
r 'k< r 'h + r 'i+ r j ■ 

We verified r' h > already. For the remainder there are two cases, according as U n = or 
U n = Sj for some i 7^ h. If U n = the three remaining inequalities in (|6.1U|) follow from 
inductive hypothesis (ii) for U'. If U n = Sf the inductive hypothesis (iii) gives the last two 
inequalities in (|6.1Uj) . and also that r\ < 0. By hypothesis (i) for U' this implies r'j,r' h ,r' k > 
hence r[ < < r' h + r'- + r' k which verifies (|6.1U|) in this case, and finishes Case 2. 
The induction is complete. ■ 
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6.2. Apollonian super-packings 

We now define an Apollonian super-packing to be a set of (ordered, oriented) Descartes con- 
figurations forming an orbit of the super- Apollonian group A s , acting on a single such config- 
uration. 

The study of Apollonian super-packings will form the subject matter of part II. We call the 
set of circles in all the Descartes configurations in such a super-packing a geometric Apollonian 
super-packing. These circles comprise four orbits of a group of Mobius transformations 

0A S ( V ) = (Sl,S 2 ,S3,S4,S] L ,S^,S^,6sf) 

where Si and were defined in §3.3. In terms of the space Aijj of Descartes configurations 
this group action is given by Qj^siV) C Aut{Qw) defined by 

Q A s{V) -W^/Wp, 

using the isomorphism Aut(Qw) = M6b(2) x {I, —1} given in Theorem 17.21 of Appendix A. 

In part II we will show that the individual circles in a geometric super-packing form a 
"packing" in the weak sense that no two circles cross each other, although circles can be 
nested. This is a remarkable geometric fact, because we will also show in part II that a strongly 
integral Apollonian super-packing necessarily contains a copy of every integral Apollonian circle 
packing, with all these copies (essentially) contained inside the square 0<x<2,0<y<2. 
All these Apollonian packings fit together in such a way that all the circles in these packings 
manage not to cross each other. 
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7. Appendix A. Mobius Group Action 



The (general) Mobius group M6b(2) is the group of Mobius transformations, allowing reflections 
(including complex conjugation). This group is denoted GM(R 2 ) in Beardon L 3, p. 23], and 
is also known as the conformal group. The group M6b(2) + ~ P5L(2,C) = SL(2, C)/{±I} 
consists of the orientation-preserving maps of C = M 2 U {oo}. The group M6b(2) is a six- 
dimensional real Lie group which has two connected components corresponding to orientation- 
preserving 5 and orientation-reversing transformations. It can be written as a semi-direct 
product 

M6b(2) M6b(2)+ x {l,c}, 

in which c denotes complex conjugation. Mobius transformations take circles to circles (or 
straight lines) and preserve angles. Thus they take ordered Descartes configurations to ordered 
Descartes configurations. Concerning orientation, we show that although Mobius transforma- 
tions can reverse orientation of a single circle, they preserve (total) orientation of oriented 
Descartes configurations. Recall that the total orientation of an oriented Descartes configura- 
tion is the sign of the sum of the (signed) curvatures of the circles in it. 

Theorem 7.1. Mobius transformations preserve the total orientation of oriented Descartes 
configurations. 

Proof. Positively oriented Descartes configurations are characterized by the four circles hav- 
ing disjoint "interiors", specified by the positive orientation. Mobius transformations take 
Descartes configurations to Descartes configurations, and preserve the "disjoint interior" prop- 
erty, hence preserve positive orientation, as given by normal vectors to the circle. The result for 
negatively oriented Descartes configurations holds since all normal vectors are reversed from 
the positively oriented case. ■ 

We now consider a group GM*{2) with four connected components, which we will term the 
extended General Mobius group, defined by GM*(2) := M6b(2) x {I, —I}. Here {I, —1} are in 
the center of this group, and we write elements of GM*{2) as ±g, in which q £ M6b(2), and 
the sign indicates which of ±1 occurs. We have 

GM*{2) s 51,(2, C) x {I,c}. 

The main object of this Appendix is to define an action of GM*(2) on the right on the 
parameter space given in the next theorem. This amounts to finding an explicit isomor- 
phism between GM*{2) and Aut(Qw), which appears as the horizontal arrow on the left in 
the top row in Table ^ This map when restricted to the smaller groups M6b(2) and M6b(2) + 
give the other two horizontal isomorphisms on the left side of the table. Table A-l also indi- 
cates isomorphisms on its right side to the Lorentz group 0(3, 1) and corresponding subgroups, 
which we defer discussing until after the following result. 

Theorem 7.2. Let GM*(2) := Mob(2)x{I, -I}. There is a unique isomorphism it : GM*(2) - 
Aut(Qw), with image elements V± s := 7r(±g) ; such that the following hold. 

5 The terminology "orientation" here refers to the invariant det(M). 
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GM*(2) = Mdb{2) x {-1,1} — ► Aut(Q w ) 



0(3,1) 



7T 



7T 



7T 



Mdb(2) = PSL(2,C) xi U,c} — Aut(QwO 



0(3,1) 



Mdb{2) + = PSL(2,C) ^ Aut(Q w )l 0(3, l) 1 



Lorentz group 



orthochronous 
Lorentz group 



proper orthochronous 
Lorentz group 



Table 1: Group Isomorphisms 



(i) For g £ Mdb(2) the augmented curvature- center coordinates for each ordered, oriented 
Descartes configuration V satisfy 



w g{v) = w v v-\ 

(ii) The action of —I on augmented curvature- center coordinates is 

= WpVlj = -W c . 



(7.1) 



(7.2) 



Proof. We compute the action of M6b(2) acting on augmented curvature-center coordinates. 

Let (b, b, w±, W2) = C 1+V r T i r' ^r) ' 3e * ne augmented curvature-center coordinates of the 
circle 

(x - xif + (y- yif = r 2 . 



This circle can be recovered from these coordinates via 

(bx - wi) 2 + (by - w 2 ) 2 = 1, 



(7.3) 



and the orientation of the circle (inside versus outside) is determined by the sign of b. An 
oriented "circle at infinity" is a line given by 



x cos 9 + y sin 9 = m, 
and its associated curvature-center coordinates are 

(b, b, wi, W2) = (2m, 0, cos 9, sin 9). 



(7.4) 



(7.5) 



Here the orientation is given by the convention that the normal (cos 9, sin 9) points inward. 
The group M6b(2) is generated by 

(1) translations t ZQ = z + zq; 

(2) dilations d x (z) = \z with A G C, A / 0; 

(3) the conjugation c(z) = z; 
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(4) the inversion in the unit circle )(z) = •= = t-tj. 

Given g G M6b(2), we will let g denote the corresponding action on the curvature-center 
coordinates of an oriented circle. The action of translation by zq = xo + iyo is 

i zo (b, b,wt, w 2 ) = (b + 2^1x0 + 2w 2 yo + b(xl + yl),b,wi + bx , w 2 + by ). (7.6) 

The action of a dilation with A = re ld , (r > 0) is given by 

T)\(b, b, w\, w 2 ) = (rb, -, w\ cos 6 — w 2 sin8, w\ sva.9 + w 2 cos 6). (7-7) 
r 

The action of complex conjugation is 

c(b, b, wi, w 2 ) = (6, b,u>i,-w 2 ). (7.8) 

The action of inversion is 

j(6, b,wi,w 2 ) = (b,b,wi,w 2 ). 
All of these actions apply to "circles at infinity" and extend to linear maps on the 4x4 matrices 

The translation operation is given by right multiplication by the matrix 



1 

x l + vl 1 x o yo 

2x 10 

2y 1 



(7.9) 



and one verifies (|7.1j) holds by direct computation. 

For the dilation with A = re 1 ® (r > 0) the right action is by the matrix 



V 



r 

l/r 

cos 6 

- sin I 





sin / 

cos ( 



(7.10) 



For complex conjugation c, the right action is by the matrix 

V7 1 = V r := 



10 

1 

1 









-1 



(7.11) 



For the inversion jc in the unit circle, the permutation matrix 



vr 



V 



Pi 



(12) 



10 
10 
10 
1 



(7.12) 



It is easy to verify that the above matrices are all in Aut(Qw)' , so that the map so far 



defines a homomorphism of M6b(2) into Aut(Q 



w 



(9(3, 1)1, identified with the isochronous 
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Lorentz group. The group M6b(2) acts simply transitively on ordered Descartes configurations, 
as observed by Wilker [I2J Theorem 3, p. 394], and the group Aut(Qw) acts simply transitively 
on ordered, oriented Descartes configurations by Theorem 13.31 Because Aut(Qw) is of index 
2 in Aut{Qw) — 0(3, 1), we conclude that the map so far defines an isomorphism of M6b(2) 
onto Aut(QwY ■ 

To complete the proof, we define the action of —I to be 

(V-i)- 1 = V_i = -I. (7.13) 

It has the effect of reversing (total) orientation of the Descartes configuration, and does not 
correspond to a conformal transformation. Since —I ^ Aut(Qw) , adding it gives the desired 
isomorphism of GM(2) onto Aut(Qw)- ■ 



In terms of the natural coordinates on M6b(2) + ~ PSL(2, C) the homomorphism ir(-) 
given in Theorem 17.21 is a non-linear map. This can be clearly seen in (|7.1U|) . where both U 0a 
and U" 1 are nonlinear functions of the coordinates of 



Da 



yf\ 



E 5L(2,C). 



(7.14) 



Here the two choices of ±y/\ give a well defined matrix modulo ±1. 

We now return to the data in Table Q giving the isomorphisms of Aut(Qw) and its 
subgroups to the Lorentz group 0(3, 1) and its two subgroups 0(3, 1)^ the orthochronous 

T 

Lorentz group, and 0(3, the proper orthochronous Lorentz group, using the terminology of 
Wilker An orthochronous Lorentz transformation is one that preserves the arrow of time; 
Table Q shows this corresponds to preserving total orientation of a Descartes configuration. 
The set of isomorphisms given by the three horizontal arrows on the right in Table ^ are 
obtained by any fixed choice of matrix A that intertwines Qw and Qc by Qiy = A^Q^A, in 
which case the isomorphism is Aut(Qw) = A _1 0(3, 1)A sending V ^ AVA -1 . To preserve 
the underlying rational structure of these forms, one must choose A to be a rational matrix, 
for example the integral matrix given in (|3.17j) . 

There is a different choice of intertwining matrix, for the rightmost horizontal rows in Table 
A-l, which will allow us to make the composed horizontal maps in Tabled from the Mobius 
group to the Lorentz groups compatible with the framework used by Wilker |42j . In |421 
Theorem 10] Wilker gave an explicit isomorphism of PSL(2, C) onto the proper orthochronous 
Lorentz group 0(3, 1)I_ It maps 6 



a 

c 



b 
d 



+ \b\ 2 %ac + bd) 



+ M 2 + M 2 ) 



+ |c| 



-ab 

,|2 



— cd) $l(ad 
-\b\ 2 9(ac- 
\d\ 2 ) 



■be) 

bd) 



5R(a6 + cd) 3f(ad + be) 



+ \b\ 2 ^{ac + bd) 
- \d?) 

Q(—ab + cd) Q(—ad + bc) 
5R(ac - bd) 

3?(aJ + be) 



(\a\ 2 -\b\ 2 
\c\ 2 + \d\ 2 ) 
3?(a6 - cd) 



We have permuted the first and last row and column of Wilker's result because he uses the Lorentz form 

2 



2 i 2 i 2 

w +x +y 
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One can identify this map with the composed map in the last row of Table A-l from the first 
column to the last column, if one chooses the map Aut(Qw) — ► Aut(Qc) to be V ^ ZVZ -1 
for a particular intertwining matrix Z, which is unique up to multiplication by ±1. It has 
irrational entries, and is 



±y/2 



110 
-1 
1-10 
0-10 



(7.15) 
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8. Appendix B. Groups in Hyperbolic 3-Space associated to 
Apollonian Packings 

In this appendix we describe two different discrete groups of isometries of hyperbolic 3-space 
H 3 associated to an Apollonian packing. We specify these groups as groups of Mobius trans- 
formations of C. 

Any Mobius transformation of C = K 2 has a unique Poincare lift to an isometry of hyper- 
bolic 3-space H 3 , viewed as the upper half-space of M 3 = M 3 U {oo}, as described in Beardon [3J 
Sect. 3.3]. Here C is identified with the ideal boundary of M 3 . In this way any group of Mobius 
transformations acting on C lifts to a group of hyperbolic isometries. 

A classical Schottky group is a group r = (fli, ...Q n ) of Mobius transformations associated 
to In circles in the plane C having disjoint interiors, but possibly touching on their boundaries. 
The group has n generators Qj, in which the j-th generator of the group maps the exterior of 
circle C%j—\ onto the interior of circle C%j. Schottky groups consist entirely of holomorphic 
Mobius transformations. (Note: knowing a single Qj alone does not determine the circles C<ij~\ 
and C2j uniquely.) 

Mumford, Series and Wright 33, Chapter 7] observe that given an Apollonian circle packing 
Vt> generated by a Descartes configuration T>, there is a Schottky group T(T>) in hyperbolic 
3-space H 3 whose limit set is exactly the limit set A© of the Apollonian packing Vt>- The 
Schottky group T(T>) has two generators Pi,P2, involving maps between two pairs of circles 
in the dual Descartes configuration V- 1 , specified by the ordering as A,a,B,b, say. It takes 
the exterior of circle A to the interior of circle a and the exterior of circle B to the interior of 
circle b. 

We treat the case they picture Figure 7.3]. where the Schottky limit set is A(X>o), for 
the (ordered) Descartes configuration Dq having coordinates 



1 -1 
2 

2 

1 3 





1 
-1 

-2 



This configuration is pictured in Figure 1101 with part (a) indicating curvatures and part (b) 
shows the ordering of circles, used later in labeling Mobius group generators Sj. 





Figure 10: Special Descartes configuration Dq 
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The Schottky group operations are associated to the dual Descartes configuration T>q, 
which has augmented curvature-center coordinates 





1 1 
1 1 

4 




-1 
1 




It is pictured in Figure ITT1 with (a) indicating curvatures and (b) the labeling of circles used 
in the Schottky group. 









Figure 11: Dual Descartes configuration Vq 



In Figure ITT1 the Schottky group generator Pi maps the exterior of circle A to the interior 
of circle a, and P2 maps the exterior of circle B to the interior of circle b. We let Ts(T>) 
denote the Schottky group they generate. The limit set Ar s (x>) of this Schottky group is the 
closed set in P 1 (C) obtained by removing the interiors of all circles in the Apollonian packing 
associated to V. This is exactly the residual set of the Apollonian packing, see Theorem 14.21 
for its properties. It is a fractal set of Hausdorff dimension about 1.3. 

The precise result is as follows. 

Proposition 8.1. The Schottky group V = (PijPg) given by the parabolic generators 



l-i 1 
1 l + i 



and 



1 

-2i 1 



U) 



has a limit set Ar which is the complement of the interiors of the circles in the Apollonian 
packing Vt> associated to the Descartes configuration T>o above. 

This limit set, "the glowing gasket", is illustrated in Figure 7.3 in Mumford, Series and 
Wright [23 Chapter 7]. 

The Schottky group Ts(T>q) is very special; it admits no deformations other than conjugacy, 
and under conjugacy 

grs(£> o )0 -1 = r s (£>) with V = q(V ). Recall that a matrix M G 
PSL(2, C) is parabolic if Trace(M.) = ±2; its action on the Riemann sphere has a single fixed 
point. Here Pi,P2 are parabolic, as is 



[Pi,P 5 



P 1 P 2 P 1 ^ 1 P2 1 



-1 - 2i 
-2i 



2i 
■1 + 2? 



J.2) 
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The Schottky group rs(Do) is characterized up to conjugacy by the property that it is gen- 
erated by two elements Pi,P2 which are parabolic, such that their commutator [P^Pg] = 
P1P2P1 X PJ 1 is also parabolic ( 33, pp. 207-208]). 

Now consider a the group G_a{Vq) = (51,52,53,54) of Mobius transformations corresponding 
to the Apollonian group for the Descartes configuration T>q, see ([4.2)1 . Here Sj is inversion in 
the i-th circle of the dual Descartes configuration T> ± , which is the circle passing through the 
three intersection points in V not touching Cj. The group G^(Vq) lifts to a group r^(2? ) 
of isometries of H 3 . We describe Ga(T>q) for the case that T>q is the Descartes configuration 
in Figure ITU1 with circles numbered as in Figure E3(b) ; all other G^(T>) are related to it by 
conjugacy in the Mobius group. An inversion in a circle C of radius r centered at zq is the 
anti-holomorphic map 

r 2 z z + r 2 - \z \ 2 
)c{z) = zo + - — = = = — z • 

Z - Zq Z- Zq 

We can write the group generators as Sj = pj o c, where pj G P£X(2,C) is holomorphic and c 
denotes complex conjugation, with 





1 — i —i 




' 1 + i -i 




1 " 




"10" 


Pi = 


i 1 + i 


, P2 = 


i 1 — i 


, P3 = 


4i 1 


! P4 = 


1 



Since p4 = I, the subgroup of holomorphic elements of Gj^{T>) is 

G\{V) = (pip 2 , P1P3, P1P4, P2P3, P2P4, P3P4, P2P1, P3P1, P4P1 , P3P2, P4P2, P4P3) 

= (P1,P2,P3>, (8.3) 

where we used the fact that pj = pj . 

There are clearly relations between the Schottky group Ts(T>q) and the group G^(Do) 
of holomorphic elements of Ga(T>) Indeed the group G^(Po) is generated by three parabolic 
elements pi,p2,p3, and these are related to the three parabolic elements P\ , P2 , [-Pi , P2] in 
Ts(2?) by the relations 

p 3 = (P 2 )- 2 , p-% = -(Pi)" 2 , p 2 = -[P 1 ,P 2 }. 

However these groups are not exactly the same. It remains to determine a more precise relation 
between them, explaining how they are give rise to the same Apollonian packing. Note that 
both these groups lifted to discrete isometry groups of H 3 have fundamental domains in H 3 of 
infinite volume. 
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